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Q ■ Abstract 

CN ' In the very early Universe matter can be described as a conformal invariant ultra- relativistic 

perfect Huid, which does not contribute, on classical level, to the evolution of the isotropic 

^ I and homogeneous metric. However, in this situation the vacuum effects of quantum matter 

^ ! fields become important. The vacuum effective action depends, essentially, on the particle 

content of the underlying gauge model. If we suppose that there is some desert in the 

(^ ■ particle spectrum, just below the Planck mass, then the effect of conformal trace anomaly 

^ ■ is dominating at the corresponding energies. With some additional constraints on the gauge 

_^ ' model (which favor extended or supersymmetric versions of the Standard Model rather 

than the minimal one), one arrives at the stable intiation. In this article we report about 

D ! the calculation of the gravitational waves in this model. The result for the perturbation 



spectrum is close to the one for the conventional inflaton model, and is in agreement with 
the existing Cobe data. 
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1 Introduction 

Inflation is considered today as a necessary component of the cosmological standard model. 
Its actual realization can be made through numerous inflaton models (see, for example, 
||l[). Besides to present a simple solution to the flatness and horizon problems, such issues 
as metric and density perturbations have been successfully studied in the context of the 
inflationary models and led to a consistent scenario for the structure formation and for the 
anisotropics in the relic radiation. These theoretical studies produced sufficient amount of 
information that could be checked in the observations and experiments, including the recent 
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(and especially future) Cobe data. At the same time, from our point of view, there is a lack of 
a natural model for inflation. In particular, the inflaton potentials which are used to generate 
the successful inflation are, in the most of the models, postulated in some appropriate way. 
In other words, those are phenomenological potentials, which can be hardly derived from 
some quantum field theory. On the other hand, the inflaton itself should be some scalar 
field with the VEV of the Planck order. The existence of such field is inconsistent with the 
modern high energy physics, which favors other candidates, like the superstring, for the role 
of the fundamental theory. Some criticism could be attributed also to the string inflationary 
models, since they require a very special initial conditions 0. 

Since the universe expands, during inflation, for many orders of magnitude, the typical 
energy scale greatly decreases, and it is reasonable to look for a model of inflation which 
should be robust with respect to this change. An alternative approach to inflation, which 
satisfies this condition, can be based on the effective action of gravity resulting from the 
quantum effects of matter fields on the classical gravitational background P, H, §, P, 0. The 
general physical input of this approach is the following [^ |[. Suppose there is a desert 
in the spectrum of particles which extends to some orders of magnitude below the Planck 
scale. Then at these energies the adequate microscopic description is that of the effective 
low-energy quantum field theory. This theory can be the Standard Model (SM), or some 
GUT. Due to the existence of the desert, in the very early Universe matter may be described 
by the free radiation, that is, microscopically, by the set of massless fields with negligible 
interactions between them. Because of the conformal invariance, these fields decouple from 
the conformal factor of the metric. We suppose, for a while, that the metric is isotropic and 
homogeneous. In this situation, the dominating quantum effect is the trace anomaly coming 
from the renormalization of the vacuum action. The anomaly-induced effective action can 
be found explicitly [^ |T0| with accuracy to an arbitrary conformal functional which vanishes 



for the special case of the conformally flat metric ||TI| . 

Treating the anomaly- induced action [^ |I^ as quantum correction to the Einstein-Hilbert 
term, one can explore the possibility to have inflationary solutions, investigate their depen- 
dence on the initial data and perform the analysis of density and metric perturbations. In 
this letter we shall concentrate on the later and derive the spectrum of the gravitational 
waves on the background of the anomaly-induced inflation. 

The article is organized as follows. In the next section we present a brief review of the 
effective action induced by anomaly. In section 3 we discuss the shape of the inflationary 
background and establish the restrictions on the particle content of the gauge model which 
produces the vacuum quantum effects. The new aspect, as compared to the previous papers 
1^ 1^ is that we obtain the inflationary solution on the basis of covariant and local version 
of the effective action. The possible solution of the grace exit problem is also discussed. 
In section 4 we derive the equation for the metric perturbations. The restrictions on the 
initial data for the auxiliary fields and metric (vacuum state for the metric perturbations) 



are implemented using earlier results [|1^ for the black hole vacuum in the same anomaly- 
induced theory. In section 5, the results of numerical analysis of the spectrum are exposed 
and analyzed. In the last section we draw our conclusions and present discussions including 
possible future steps in investigating the model. 



2 Effective action induced by anomaly 

Since we shall need many particular details of the anomaly-induces action, it is reasonable to 
present its derivation, also in some details. The starting point is the action of free massless 
fields in curved space-time: A'"o real scalars, N112 Dirac spinors, and A^^i vectors. All A^'s 
indicate a number of fields, not multiplets. The conformal versions of the actions, for each 
of these fields are: 

S,= jd'x./^g{-\F,,Fn. (1) 

The only divergences which one meets for the free fields are the one-loop vacuum ones. 
Using the well-known results for these divergences and the notations 



for the square of the Weyl tensor, and 

E = R^vapR^ — 4 RapR + R 
for the integrand of the Gauss-Bonnet topological term, we get 
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j d'^x^^{wC'^-hE + cUR}, (2) 



where e = (Att) {D — 4) is the parameter of dimensional regularization. 
The renormalized one-loop effective action has the form 



r = ^ + r + A^, (3) 



where T is the quantum correction to the classical action. AS is a local counterterm which 
is called to cancel the pole in (j^) . The classical action of the renormalizable theory is given 
by the sum S = Smatter + S vacuum, whcrc Syacuum must iuclude the following structures: 

Svacuum = J d^x^^ [aiC^ + tts^ + agDi?} + ... . (4) 

Here 01,2,3 are the parameters of this vacuum action. The renormalization of this parameters 
is necessary element of the consistent quantum theory. If we do not include these terms 
into the classical action, they will arise anyway due to the quantum corrections and with 
unremovable divergent coefficients 0]0. The one-loop /5-functions for the parameters 01,2,3 



are given by (g). 

The vacuum action may also include some non-conformal terms like the Einstein-Hilbert 
one, cosmological term or -y/^i^^-term, but their renormalization is not necessary in the 
case of conformal invariant free massless fields which we are dealing with. But they may 
be, indeed, important from other points of view. In particular, later we shall include the 
Einstein-Hilbert action into Svacuum- 

AS in (^ is an infinite local counterterm which is called to cancel the divergent part 
(H) of r. Indeed AS* is the only source of the noninvariance of the effective action, since 
naive (but divergent) contributions of quantum matter fields are conformal. The anomalous 
energy momentum tensor trace is IT^ ITS 



with the same coefficients w, b, c as in (|). The eq. (|) can be also considered as the equation 
for the finite part of the 1-loop correction to the effective action. The solution of this equation 
is straightforward ^, |l^, 11 1. Let us consider the metric in the form g^^, = Qf^uc'^'^. Then 



the effective action is given by the expression [^, 0, [Tl[| : 

1 

(4^ 



Sclg^u] + T^ Jd^x./^{waC^ - ba{E - ^Di?) - 2baA,a - 



Here 



-^ic-h)[R-6{Vaf-6ina)r)}. (6) 

A4 = □' + 2 R'^'V^V, - ^ i?n + 1 (V^i?) V^ 

is the fourth derivative, conformal invariant and self-adjoint operator and Sc[gfj,u] is some 
unknown functional of the metric ^^^.(x), which serves as an integration constant for any so- 
lution of (^). The action (^ includes some arbitrariness, which were extensively investigated 

^The importance of these high derivative terms concerns also any other inflationary model which is going 
to deal with the quantum fields. 
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recently [0, |T4|. Of course, all the arbitrariness is inside the conformal functional Sc[gi_iv\- 
If one succeeds to rewrite (^ in terms of the original variable g^^, this functional should be 
substituted by a conformal-invariant functional of this metric. 

The expression (P) can serve as a basis for the inflationary solution 0. This solution 
is identical to the one found by Starobinsky [|] using the (OO)-component of the corrected 
Einstein equations. However, since we are going to consider the perturbations of the metric, 
it is better to maintain covariance. Therefore one has to rewrite (|^) in terms of the original 
metric g^^,. The action (^ can be presented in a nonlocal but covariant form using the 
original metric and then in a local covariant form via an auxiliary scalar H]. We shall 



perform this, following the scheme developed in [Tj] and applied to the derivation of the 
Hawking radiation of the black holes in |T^. The non-local form for the anomaly-induced 
action is 

-g{x)R\x)+ 



Sc[g, 



lb 



/lU] 



+ 



12(47r) 
2 



d X 



J d^x^-g{x) J dS^-g{y) {E - -UR), G{x, y) 



W^2 

4^ 



^-H 



Sc[g, 



c-ib 



flUj 



12(47r) 



d'^x^-g{x) / d^y^-g{y) 



xG{x,y) 



Vb 



2 , d'^x\J —g{x) R^{x) 



Vb 



iE-luK)-^-C^ 



(^-|ofi)-^C^ 



+ 



X 



+ \jd'x4^^) jdS^^^)(^C 



G{x,y) 



w 



2Vb 



C' 



(7) 



Here G{x, y) is a Green function for the operator A4. One has to notice that we have 
introduced an additional (as compared to 0) structure J C^ G C^ in order to write the first 
non-local term in the symmetric form. The importance of this term to be included into the 
conformal part of the effective action Sc[g^i,] has been previously discussed in [|1^, |T5|] and 



recently in the last of Ref . |T6| . 

The second term in (^ is local, but the others are not. However the last two terms can 
be done local through the introduction of the auxiliary scalar fields. Thus we arrive at the 
following final expression for the anomaly generated effective action of gravity. 
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The last action is classically equivalent to (^, for if one uses the equations for the auxiliary 
fields ^ and ■?/', the nonlocal action (0) gets restored. At the same time, local theory 
(^) is much more useful for the applications. The action (^, exactly as other forms of 
the anomaly-induced action, contains the arbitrariness related to the conformal invariant 
functional S^g^J\. One has to notice that the only relevant classical term J wC"^ in the 
classical action of vacuum (^) is conformal invariant and it can be unified with ScIquu]- This 
conformal invariant functional is not relevant when we are interested only in the behaviour 
of the conformal factor of the metric. In this case the result (|^) is exact one-loop correction 
to the effective action. At the same time, any other application of (^) requires fixing this 
arbitrariness in this or that way, so the consideration becomes, in part, phenomenological. 
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For instance, in 15 the conformal functional has been set to zero, and this led to the 



classification of the vacuum states for the semi-classical black holes and to correct, in the 
leading order, result for the Hawking radiation. Therefore, the Sc[g^u] = choice can serve 
as a reasonable approximation and in general we adopt it here. 

The action (H) contains high derivative terms and some remarks are in order. 

First: our purpose is to investigate the classical equations of motion following from the 
anomaly-induced action. That is why here we do not perform the path integration over the 
auxiliary fields (p and ip- As a result the values of the coefficients w, b, c remain unaltered, 
exactly as it was in the similar black hole application |T^, but in contrast to the original 
work [P]. 

Second: the kinetic term for the auxiliary field ip is positive while for ip it is negative. 
This indicates that these fields should not be considered as physical, but only as auxiliary 
ones. 

Third: the very fact that the above action contains higher derivatives does not mean 
that it can not be applied to the consistent description of physical phenomena. Let us 
remind that we do not consider the quantum theory of gravity, and metric here is nothing 
but the classical background 0. The higher derivatives show up only in the vacuum part and 
thus do not jeopardize the unitarity of the quantum theory 0. In our opinion, for classical 
theory the only one reasonable criterion of consistency is the existence and stability of the 
physically acceptable solutions. As we shall see in the next section, the anomaly-induced 
effective action really produces such solutions. 

As it was already mentioned above, for the quantum field theory of matter fields in the 
curved space-time the high derivative classical action of vacuum (H[) is necessary, because 



otherwise the theory is inconsistent [|Tl|. Thus, the appearance of high derivative quantum 

^One has to notice that the treatment of the corresponding quantum corrections as the 1/N approximation 
to quantum gravity |17, m pi is inconsistent Mi 



^At the same time, there are certain indications to that even the high derivative quantum theory can be 
indeed unitary |T9] . This was recently found to be true for the particular version of the four-dimensional 
anomaly- induced theory which we are investigating here pOf . 



corrections in (P) does not change this aspect of the vacuum action. On the other hand, 
the high derivative parts of the vacuum action (classical and quantum) can be important 
only in the high energy domain, while at the low energies they should be treated as a very 
weak correction to the Hilbert-Einstein action. Therefore, the cosmological application of 
the above action is essentially restricted by the short time after Big Bang, in which the 
typical energy of matter does not decrease too much. However, in this short time inflation 
happens, and we consider this in the next section. 

Before going on with the study of cosmological applications we shall comment on the 
possibility to reduce the action (^) to the action without the higher derivatives via the 
properly chosen additional auxiliary fields. In general, this can be done using the trick 



suggested in ||2T[ but, unfortunately, it is not really helpful for the cosmological applications. 



For instance, introducing two extra scalars x ^i-iid 6', one can present (^) in the form 



{At:) 



l_2 , Vb ^ 1 . , V3^^^2b 



+ Jd'xV^i --X' + j^xO^ + 2^ + 4^ ^ [aa+{Vay\^ . (9) 

The last equation resembles the actions of gravity coupled to two scalars (which could be 
even called infiatons). So, one can try to reduce the study of the cosmological consequences 
of the above action to the standard procedure of dealing with infiatons. However, this 
program meets serious difficulties. The most explicit one is that the reduction of the order 
concerns only a sector. The tensor degrees of freedom are still with fourth derivatives, and 
for their reduction one has to introduce tensor auxiliary fields. This means, that there is no 
explicit way to use the standard infiaton-based results in our case. At the same time, this 
presentation indicates to the qualitative similarity between two approaches, and gives some 
hope that the main prediction of the theory (^ will not be very different from the ones of 
the phenomenological inflaton models. 

3 Inflationary background 

In this section we shall consider the inflationary solution for the dynamical equations of the 
theory with the action 

r2 



Stotai = -M'p J d^x y^R + r , (10) 

where Mp = I/IGttG is the square of the Planck mass, and the quantum correction T is 
taken in the form (^. In what follows we set ScIQiiu] = 0, including to it, when it is not 
indicated explicitly, also the classical vacuum term. 



Since we are going to look for the isotropic and homogeneous solution, the starting point 
is to choose the metric in the form g^^ = a?{ri)g^y, where rj is conformal time. It proves 
useful to denote, as before, a = In a. Now, one has to derive the equations for the three fields: 
yj, ^, and a. In the rest of this section we shall consider the conformally flat background 
and thus set g^y = r]^^. 

The equations for ip and tp have especially simple form 
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One has to remind the transformation law for the quantities which enter the last expression: 
y/^C = V^C^ , v^ A4 = v^A4 , (11) 



-g{E--nR) 



-g{E--aR + AA^a). 



(12) 



Taking into account our choice for the fiducial metric g^^i, 
equations in flat space-time 



rjn,y, one arrives at the following 
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n2^ = 0. 



(13) 



The solutions of ([T3| ) can be presented in the form 



^ 






ijj = 1pQ. 



(14) 



where y^o, '^0 are general solutions of the homogeneous equations 0^(^0 = 0, U^ ipQ = 0. 
Thus one meets an arbitrariness related to the choice of the initial conditions for the auxiliary 
fields if, ip. But, the inflationary solution does not depend on y^o, V'o- Substituting (|T^) back 
into the action and taking variation with respect to a we arrive at the same equation for a 
that follows directly from (^. We shall write this equation in terms of the physical time t, 
defined, as usual, through a{ri)dri = dt. The useful variable is H{t) = a{t)/a{t) = o"(t), since 
the equation (completely equivalent to the one of 0) is of the third order in this variable: 



H + 7HH +41+ 
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(15) 



It is easy to identify the special solution corresponding to if = const: 

Mp 
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Vb 



a{t) = ao ■ exp Ht . 



(16) 



Positive sign corresponds to inflation. The solution (|T6D had been first discovered in [^, |^ 
and studied in |, || [|. 

The detailed analysis shows |^, 0] that the special solution (|TB]) is stable with respect 
to the variations (not necessary small) of the initial data for a{t), if the parameters of the 
underlying quantum theory satisfy the condition - > 0, that leads, according to (H), to the 
relation 

iVi < ^iVi/2 + ^iVo. (17) 

This constraint is not satisfied for the Minimal Standard Model (MSM) with A^^i = 12, 
A'^i/2 = 24 and A^o = 4. However, one can consider some aspects of the neutrino oscillations 
as an indication that the MSM should be extended, and in this case the inequality (|17D can 
be readily satisfied. Below we consider two versions, each of which leads to stable inflation. 

i) Extended SM with A^i = 12, Ni/2 = 48, A^q = 8 

and 

ii) Supersymmetric MSM with A^i = 12, N1/2 = 32, A^'o = 104. 

The advantage of stable inflation is that it occurs independent of the initial data. After 
the Big Bang, when the Universe starts to expand and the typical energy decreased below 
the Planck order, we can imagine some kind of " string phase transition" . Starting from this 
point, the effective quantum field theory is an adequate description, and the anomaly- induced 
model applies. In case of the stable inflation, the initial data for a{t) and its derivatives do 
not need to be fine tuned, if only the condition (|T7D is satisfied - the inflation is unavoidable. 

Let us calculate the necessary duration of inflation. Suppose we want the Universe to 
expand in n e-folds. Then the total rate of inflation (in the Planck units of time) will be 

^^^^^±^ = exp|4vr./|^At|, AT,,, = A^o + 11 ■ AT^/, + 62 • AT, (18) 



so that 
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For the extended and supersymmetric versions of the Minimal Standard Model the time 
necessary for 65 e-folds is around ten Planck times only. The numerical study has shown 
that the exponential solution stabilizes, in the theories satisfying ([T7|), in much shorter time. 
For that reason, one can safely derive the metric perturbations on the exponentially inflating 
background, independent on the initial data. At the same time, the numerical analysis has 
demonstrated, that the non-stable inflation favored in the original works of Starobinsky 0, |^ 
may be inconsistent from the cosmological point of view. The system leaves the inflationary 
phase too fast and the necessary 65 e-folds can not be achieved in the a ~ exp{Ht) regime. 

^In Q two other similar solutions for the FRW metric with fc = ±1 were found. 



The most difficult question is how the inflation ends. So far, we do not have a definite 
answer to this question, but there are some particular indications that a solution is possible 
if we take the masses of the matter particles (perfect fluid) into account. We can mention, 
at the ffist place, that the ffist study of the same model, with density of matter p ~ a~^ 
inserted into the (OO)-component of the equations, has demonstrated the transition to the 
FRW behaviour at the late time limit 0. This shows, at least, the possibihty of a desirable 
particular solutions due to the matter fields. 

A very important observation is that even when the Universe expands so rapidly as 
in (p!6D, the gas of matter particles performs some work, and the average energy of these 
particles decreases. At some instant this energy decreases such that their masses become 
relevant and then the matter part of the equation obtains some dust component. It is 
easy to see that in this case (|16[) is not anymore a solution. The classical solution for dust 
a{t) ~ t^/^ also is not a solution because of the quantum term. However, both Einstein and 
matter terms behave like t~^ while the "quantum" part behaves like t~^, and very rapidly 
the anomaly-induced quantum term becomes irrelevant. Thus, one can suppose that at the 
long-time limit a{t) ~ t^/^ is a good approximation for the unknown solution of the equation 
with matter. Of course, the above consideration is not a complete solution of the grace exit 
problem, but one can hope to find such a solution along this line. In the rest of this article we 
will not discuss the grace exit problem, but instead concentrate on the gravitational waves 
and their spectrum during the inflationary period. 

4 Gravitational wave equations 

The equation for the metric perturbations is based on the bilinear expansion of the action 
of interest (plQl). Before going on to this cumbersome expansion, let us present the action, 
through some integrations by parts, in a more convenient form S = J d^x L , with 

L = J2fsLs = V^\ foR + fiR'^^^'Raf^-ys + f2R''^Rai3 + fsR^ + U^aR + f,^Aip\ . (19) 
Here the following definitions have been introduced: 

/o = -Ml, ; (20) 

b — w w 

h = ai + a2 + - — 7f^ + - — if^ ; (21) 

h = -2ai-4a2 + ^^^¥^--^V' ; (22) 

ai 3c — 2b 3b — w "^ , 

10 



f^ = I ■ (25) 

One could simplify the coefficients /4 and /s, substituting their values, but it is better to 
keep them in order to track back the origin of each term in the final equations. Now we have 
to fix the arbitrariness related to the homogeneous solutions ipo and ipo in ([T^). Let us remind 
that the choice of the initial data for ip and ip defines the vacuum state for the perturbations. 
One can make a useful comparison with the vacuum of the black hole background. Within 
the semiclassical approach to the black hole radiation, the vacuum which provides a smooth 
transition to the Minkowski vacuum at the space infinity is the Boulware one. Let us suppose 
that the proper cosmological vacuum for the expanding Universe reduces to the Minkowski 
one at infinite time. For the homogeneous and isotropic solution the equation o"^ ipo = 
reduces to ipo = 0, and the solution will depend on four integration constants. But, if 
one requires the finite behaviour at infinite time, the only choice is (po = const. Now, if 
requesting the correspondence with the black hole choice at the space infinity |[T^, the only 
possibility is to set ipo = 0, and do the same with ip. Therefore, the consistent choice of the 
solutions for the auxiliary fields, which we shall use in the rest of the paper, is 

/b 
^ = -^a, i: = 0. (26) 

The calculation proceeds by taking perturbations such that 

9f^iy = 9% + h^^ , (27) 

where g^ are background inflationary solutions (|T6|) 

(7j, = (l, -5,,e^*), ^ = ^' /^ = 0,1,2,3 and ^ = 1,2,3, (28) 

and h^y are the perturbations around them. Since we are interested in the gravitational 
waves, we can retain just the traceless and transverse part of h^j^, which are the purely 
tensor modes 0. Hence, the metric perturbations are submitted to the restrictions: 

d^ h'^ = , hkk = 0, (29) 

Besides, the synchronous coordinate condition, h^Q = 0, is imposed. 

The details of the bilinear expansion are postponed to the Appendix. In fact, most of 
the expansions are identical to the ones performed by Gasperini in [^| , who has studied the 



metric perturbations for the string- induced action. We have just checked these expansions 
(and found them absolutely correct). Some other terms are typical of the induced model 

^Scalar perturbations in the anomaly-induced model have been first studied in ||22| and, in the modern 
framework, in p3[. 
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under consideration and we expanded them for the first time. The final result for the 
expansions of all structures Lj/j looks like 



a^fol 3H^h^ + hh + AHhh + h^ - ^'^\+0{h^^ 

4 4 a^ 



(30) 



Li = a^fi\2H'^h^-4:H^hh-6H^h'^-l6H^hh + h^ + 4:Hhh + —\/'^hV'^h + 

(Jj 



2h^ + (H'h - 2Hh)^}+0{h 



(31) 



a'f'. 
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4 2 4 
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(33) 
(34) 
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l^'hh - In^^hh - In'^'h' - ^^02 \+0{h'' 
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(35) 



Disregarding the higher order terms, the equation of motion can be calculated as the La- 
grange equation for the action (p!9[): 

rf2 dL 



, ^2 dL _ d^2 dL 

df^ dh dV^h dt dW^h dt Oh ' dVh ' dh 



ddL ^ dL dL ^ 

Vt^^^ + tt^ = 



(36) 



The resulting equation for the metric perturbations on the background of the exponential 



inflation ([T6| ) has the form 

= /i 
+ h 



^H^'h + \H^^^h + 3i/Vo - i^'(8/i + 12/2 + 48/3 - 9^/4) 
i^'(-9/2 - 36/3) + ^-'PH^U - ^H^'h + 
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(37) 



At this moment, one has to fix the functions /s(t) as it was discussed above. Then the eq. 
provides the appropriate basis for the analysis of the metric perturbations. 



5 The initial data problem and spectrum analysis 

Even for the special choice of /s(t), the equation (^) remain quite difficult to solve analyt- 
ically, but it admits an efficient numerical study. A numerical analysis requires some care, 
mainly in two aspects: the equations must be dimensionless and the initial conditions must 
be consistently chosen. For the first point, there is no problem in the equation written above: 
if one sets the Planck mass equal to one, time is automatically measured in the Planck units. 
For the initial conditions, we consider that the perturbations have a quantum origin: the 
seed of the perturbations are quantum fiuctuations of the primordial fields. This is fixed in 
the following way. The equation for gravitational waves reduces to the equation for a scalar 
field - the coefficient in front of the tensor mode. Then the perturbations, which originate 
from the fiuctuations of the zero point energy of the quantum fields, have the spectrum 
characteristic of a scalar quantum field in Minkowski space. This "vacuum state" is well 
known []T3| : 



h{x,r]) = h{r]) e 



±m.x 



h(r]) oc 



\/2n 



(38) 



In these expressions we employed the conformal time, since with it the FRW metric becomes 
conformal to the Minkowski metric in fiat space; n is the wavenumber vector. After fixing 
this initial spectrum, one can derive how the initial amplitude depends on n. In our case, it 

becomes 

In 

2 



ho oc 



1 



he oc 



ho oc 



n 



3/2 



ho oc 



n 



5/2 



(39) 



/2n ' " V 2 ' " V2 V2 ' 

In the last expression, the derivatives are taken with respect to the cosmic time, while the 
initial conditions are given in terms of the conformal time. Taking oq = 1, the transforma- 
tion of the initial conditions into the physical time just introduces the constants of the order 
of one. After all, the form of initial conditions tells us how they depend on the wavenumber 
value. 

Using the dimensionless equation and taking the initial conditions according to the initial 
quantum spectrum (^), one can integrate the fourth order equation (|37| ) numerically. It 
is not difficult to plot the final spectrum, but first one has to decide which quantity is 
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interesting to evaluate. A crucial number is the power spectrum of the perturbations, which 
is defined as follows. Let us suppose that the desirable solution is given by the function 
h{t,x). Through a Fourier transformation we get 

hn{t) = ^^Jhit,^)e^'^-d'x . (40) 

The square of the total amplitude of the gravitational perturbations is obtained through 

h\t) = j hl{t)£n , (41) 

and this can be rewritten as 

h\t) = 47r f hl{t)n'^dn = An f hl{t)n^d\nn = f P^{t)d\nn . (42) 

The quantity -P^(t) = h'^{t)n^ is called the (square) power spectrum and tells us how the 
amplitude of the gravitational waves varies in the interval between Inn and ln(?T, + dn). 

Thus, in order to calculate the power spectrum we must square the value for the gravi- 
tational perturbation at a given moment of time and for a fixed wavenumber, then take its 
logarithm and see how it varies with the logarithm of the wavenumber itself. Practically, 
after performing numerical integration, we have to plot 

lnn^/i^(t) X Inn. 

The power spectrum essentially says how the amplitude of the perturbations depends on 
their wavelength. A fiat spectrum, the Harrison-Zeldovich one, establishes a "democracy 
principle" : all perturbations reenter into the Hubble horizon during the radiative or matter 
dominated phase (after have left it during the infiationary phase |^) with the same ampli- 
tude independently of their wavelength. One has to notice that, in general, we obtain an 
expression for the small variations of n in the long wavelength limit (which are the most 
important for cosmology), when the dependence can be taken as P^ o: n^. This distribu- 
tion tells how the amplitude of the perturbations depends on n. The coefficient k is called 
spectral indice of the perturbation. 

In order to compare our results with the traditional ones, consider infiationary scenario 
based on the Einstein's equations with a cosmological constant. This cosmological constant 
can arise from some infiaton potential. Then, one meets the deSitter Universe, and the 
perturbations behave as |2^ 



Kv) = Vvc±{n)J^-^{nri) , (43) 

where c± are integration constants which may depend on n. This dependence must be fixed 
by the initial spectrum. But, if choosing the vacuum state as described above, we find that 
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those constants do not depend of n. Using the long wavelength limit approximation n -^ Q 
and taking the dominant mode in the above expression, we find 

Pn oc n^/^"^/^ = constant . (44) 

Hence, the traditional infiationary scenario predicts a fiat spectrum, with A; = 0. It is most 
relevant, that the analysis of Cobe, Boomerang and Maxima observational programs favor 
the fiat spectrum too fg^, |2|]. 

One can use these results to gauge our numerical procedure. Fixing the initial spectrum 
according to (|39D , integrating the equation for the gravitational wave for the traditional 
infiationary scenario, 

h + \— -2-\h = Q , 45 

a Va^ a] 

with a{t) = e^*, and using the numerical procedure described above, one meets, for the (^3[) 

case, 

fc^O.Ol. (46) 

It is easy to see, that this numerical result is close to the analytical one, which is zero. We 
have considered a variation of n between zero and one. This leads to initial perturbations 
whose scales are of the order of the Planck length. Applying now the same procedure for 
our model, with ai = a2 = and any ipo < 1, and considering the multiplet of the extended 
Standard Model, we found 

fc~-0.01, (47) 

which is qualitatively in agreement with a fiat spectrum. It is important to notice, that this 
result depends (although not dramatically) on the number of e-folds during the infiationary 
phase. Here, we have fixed that the infiationary phase lasts ten Planck times, leading to 
approximately 65 e-folds expansion. These were the same values employed in testing the 
numerical procedure in the traditional infiationary case exposed above. 



Some recent analysis indicates to —0.15 < k < 0.16 ^j. Hence, our model predicts a 
spectral indice different, but not very far from the one of the traditional infiationary scenario. 
Besides, the spectral indice in the anomaly-induced infiation has quite a good agreement with 
the observational data. 

Some additional consideration is in order. Our result is essentially based on the number 
of suppositions about the vacuum state. In particular, we have neglected the conformal 
invariant functional Sc[g^v]- This may be justified by the obvious success of the similar 



procedure in the black hole case |]15|. At the same time, it is interesting to check whether 
the dependence on this supposition is really strong. The most natural is to consider the 
nonzero coefficient ai in the classical vacuum action (|^) since, from the quantum field the- 
ory point of view, it is not possible to set ai exactly zero. The reason is that ai is subject of 
renormalization and the consequent renormalization group running (see, for example, [pJ]|). 
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We have performed the numerical analysis and found that for ai 7^ 0, different results for 
the spectrum may be obtained. For example, with ai < 1 and all other entries the same as 
before, we find the following relation between the spectral indice and the value of ai: 











Table 1 










ai 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


k 


-0.01 


-0.1 


-0.1 


-0.1 


-0.2 


-0.2 


-0.2 


-0.2 


-0.2 


-0.2 



As one can see from this table, when the value of ai increases, the spectrum becomes more 
and more negative. It means that, taking into account ai, the amplitude of gravitational 
perturbations decreases as the scale increases. Taking into account the numerical value of 
the /3-function w/iATrf ^ 0.02 for ai in (|) for the extended SM and for the MSSM, we 
conclude that the admissible value ai < 0.4 holds under quantum corrections. From the 
above table, it is possible to see that this value is within the observational limits. 

Among other factors, the dependence on the multiplet composition does not seem to be 
very important. At least, the spectrum is almost the same when one takes the multiplet of 
the extended Standard Model or the Minimal Supersymmetric Model. 

On the other hand, the flat (or almost flat) spectrum found for the Oi = essentially 
depends on the fact that our stable version of the anomaly-induced inflation develops only 
de Sitter like behaviour. Our model does not admit the power-low inflationary solutions 
(if it would be possible, the perturbed equations should also be modifled with the inclusion 
of terms containing H, H and so on). As we have already pointed out in section 3, the 
stabilization of the exponential inflationary background performs very fast , and one can 
suppose that during the stabilization period the initial quantum perturbations have no time 
to change their spectrum. One can consider this period as a small modiflcation of the initial 
data for the perturbations. On the other hand, the transition to the FRW solution may 
include the period of fast but power-like expansion. At the same time, since we have no 
mechanism to desribe such a process yet, we can suppose that this transition also has been 
performed rapidly, so that the perturbations do not suffer too big modiflcations. In principle, 
since the two scales Mp (corresponding to inflation) and the energy density of matter at later 
epoch (which deflnes the grace exit to FRW) are very different, it is natural to expect a fast 
transition which does not produce serious modiflcations on the indice k. Usually, a power-like 
inflationary behaviour leads to the decreasing spectrum A; < ||2^ and then, qualitatively, 
we have the same effect which is favored by ai 7^ 0. 

Let us now investigate the dependence on the choice of the initial data. In order to 
do this, one can substitute other sets of initial conditions, different from (IS^). The flrst 
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alternative choice was: 

ho oc 0.001 , ho (X 1 , ho (X 1 , Aooc 1 , (48) 

for any value of n. The numerical analysis shows that these initial conditions lead to the 
result: 

k = -0.00002 . (49) 

The second choice was 

/io = 0.001 , ho = -l , ho = 2 ho=0.l . (50) 

Then the result for the power spectrum is 

k 0.00002 . (51) 

Similar results have arisen for some other choices we tried. Thus, the result is not very 
sensible to the choice of the initial conditions, and one can expect the power spectrum very 
close to zero. This universality is a very good point, indeed. 

The information we have obtained concerns how the amplitude of the perturbations asso- 
ciated with the gravitational waves (tensor mode), depends on n. In the figures 1 and 2 we 
display the graphics for the multipole coefficients of the two-point correlation function of the 
cosmic microwave background radiation, measured in terms of the temperature fluctuations. 
This coefficient is given by the expression 

AT " 

— = Y.CiPiicos9) , (52) 

-^ 1=2 

where Ci are the multipole coefficients, Pi{x) is the Legendre polynomial, 6 is the angle 
between two observation points in the sky and ^ is the temperature fluctuation between 
these two points. The sum begins with / = 2 because we exclude the dipole momentum due 
to the motion of earth. 

Using the CMBFast code [BDl we plot the spectrum of anisotropy of the Cosmic Microwave 



Background due to gravitational waves for some of the cases specified above. In figure 1, we 
use k = —0.01, which is the prediction corresponding to ai = 0. In figure 2, we plot the case 
oi k = —0.2, that corresponds to ai = 1. In plotting these graphics there are two essential 
inputs: the power spectrum k and the matter content of the Universe. The power spectrum 
has been described above, and for the matter content we took the one which is the most 
accepted in the hterature: 5% of baryonic matter, 35% of cold dark matter and 60% due to 
a cosmological term. The figures show a small difference in the amplitude of the spectrum 
due exactly to the fact that the amplitude decreases slight differently in each case. 

Qualitatively, the predictions of our anomaly-induced theory agree with the traditional 
inflationary De Sitter results. In order to have a more detailed description, which could 
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definitely distinguish the predictions of our model from the ones based on the infiaton, one 
should analyze the scalar perturbations and the baryogenesis during the inflationary phase. 
However, the solution of this problem (see [^, |23| for the previous studies) is outside the 
scope of the present paper, which is devoted to the gravitational waves. 

6 Conclusions and Discussions 

At very high energies, when the inflation is supposed to occur, the effective mass of the 
particles may be considered zero, and a conformal description is possible. Then, the quantum 
effects of the matter fields generate the anomaly- induced effective action of the vacuum. In 
the preceding works 0, ||, the stable inflationary solution for this anomaly-induced model 
was set out. Contrary to the infiaton models, the anomaly-induced inflation comes from the 
very natural framework, for the degree of the phenomenological input is much smaller than 
in the infiaton models. In particular, the inflation does not require any fine-tuning of the 
initial data. From our point of view, it makes this approach very promising. 

The great advantage of the infiaton models is that they are very well elaborated, and 
come, under certain phenomenological suppositions, to the agreement with the available 
observational data. In order to reduce the gap between the two approaches, in this paper we 
have investigated the spectrum of the gravitational waves in the anomaly-induced model. 

The main difficulty of the anomaly-induced model (exactly the same concerns some other 
models [^) is the natural end of inflation. There are definite indications, however, that this 
problem can be solved in the framework of an effective approach, which is supposed to 
provide the dynamical mechanism terminating inflation. 

In this paper we mainly concentrated on the spectrum of the gravitational perturbations. 
In the traditional inflationary scenario the spectrum of tensor perturbations is flat, so that the 
amplitude of gravitational waves does not depend on their wavelength (so-called Harrison- 
Zeldovich spectrum). Strictly speaking, this result applies to the long wavelength regime 
only. In a quasi-De-Sitter space-time generated, for example, in the slow roll over models 
of inflation, the spectrum is quasi-flat. The observational results obtained until now, also 



agree with a flat or quasi-flat spectrum ^, ^, ^ ^ . 

In the anomaly-induced model, the analysis of the spectrum of the gravitational waves 
is more involved, because the equation for the perturbations is of the fourth order, and 
with coefficients depending on the particle content of the underlying theory and on the 
vacuum state. The particle contents considered here were those of the extended SM and of 
the MSSM. Even in this case, an analytical solution is not possible for the perturbations 
of the tensor mode of the metric, and the numerical analysis was performed. The initial 
conditions for this numerical analysis were suggested by a quantum mechanical mechanism 
for the generation of the perturbations. This led to the spectral indice compatible with the 



available observational data on the anisotropy of the relic radiation and very near the values 
obtained in the inflaton inflationary models. 

However, we must stress that it is generally supposed that the observational data com- 
ing from the anisotropy of cosmic microwave background radiation are due to density per- 
turbations linked with scalar mode perturbations. These density perturbations display a 
Harrison- Zeldovich (or quasi-Harrison-Zeldovich) spectrum, and the tensor perturbations 
are supposed to have the same features. In the standard inflationary scenario the spectral 
indices for the scalar and tensorial perturbations are connected by the expression kx = ks — 1 
p5| . In fact, the interpretation of these observational data are strongly model-dependent, 
and only crossing informations from different observational programs like CMB anisotropy, 
high redshift supernova and gravitational lensing, it may be possible to have trustful results 
for the observational parameters. On the other hand, it is expected that in the near fu- 
ture a direct measurement of the effects of gravitational waves in the anisotropy of cosmic 
microwave background will be possible, mainly due to polarization of the background pho- 
tons, and the contribution to CMB anisotropies due to gravitational waves will be separated 
from the contribution due to density perturbations. Then one will have more chances to 
distinguish the inflationary model which flts better with these new data. 

Many other aspects of the anomaly-induced gravity model presented here deserve further 
study. Besides the grace exit problem, the fate of density perturbations should also be 
investigated in the present framework. Another interesting aspect is the mechanism of 
reheating, which can be, hopefully, found along with the grace exit solution. The results 
obtained so far indicate, in any case, that the anomaly-induced model is a very promising 
candidate for obtaining a self-consistent inflationary phase. 
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useful conversations. I.L.Sh. thanks A. Belyaev for the discussion of the super symmetric 
SM. Authors are grateful to CNPq for the scholarship (A.M. P.) and grants (J.C.F and 
I.L.Sh.). 

Appendix 

To derive the equation for the gravitational waves, one has to consider the metric per- 
turbations in the Lagrangian and retain the bilinear part of it, so that to get linear terms 
in the the fleld equations. In order to be able to compare the technical details, we flx our 
notations equal to the ones of [^]. In the same paper one can flnd most of the necessary 
intermediate formulas, like the expansions for the components of the curvature tensor. We 
shall not present these expansions here, and do so only with the (/^-dependent terms which 
one does not meet in pl|. 
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The background solutions are such that 



9oo 



1, g,, = ~a'it)5,, = -e'^^'S,,. (Al) 

The value of H is constant, and we impose this condition from the beginning just in order to 



simplify the expressions (actually, our calculation, just as the ones of ||2^, were performed 
for an arbitrary H). It is useful to chose, as dynamical variables, the mixed component of 
the perturbation, which will be denoted as M = h. 

The tensor mode of the metric perturbations is defined by the relation 



9iiu -^ 9tiu + Sg/iu , where 



J^f i^fiu 



{A2) 



parametrized by transverse traceless tensor h^j,^. Besides, we fix the coordinates by the 
condition /i^o = 0. Expanding the contravariant components of the metric tensor, at first 
and second order in h we get 



§{l)g,^U ^ _|^^J,U 



§(2)gf.u ^ ^M"/,^ 



where ^'■"^ denotes the n-th order in the expansion of the corresponding quantity in powers 
of h. Similarly, for the determinant of the metric we get 



Using the intermediate formula 



s^''>^/^ = -\v^h^''h,,. 



1 . 



□y? =(p +3H if -^hhV+0 [h^ 



we arrive at the following expansions 



<5(2)(v/^(n^y 



-h^[^+3H^]^ 



hh(^(p +3H ^) ^ 



(AS) 



(AA) 



iA5) 



S(2) 



-gR^'V^^V,^ 



a^ ^ 



4 2 4 



(A6) 



5(2) 



-(^i^V^*^^^^ = c? (p" \'iE''h^^ hh + ^h +4:H hh- -h^ 
J 4 4 a"' 



(A7) 



These expansions have been used, together with the ones of ||2^, in the main text of the 
paper, for the derivation of expressions (P5|). 

Note added. After this peper has been submitted for publication, we have learned, from 
the preprint |]3^, about the very important paper [^, where the Starobinsky model has 
been analized from the point of view different from the one of H, 0. In particular, it is 
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pointed out that the necessary duration of inflation in the non-stable case can be achieved 
only through the fine-tuning of the coefficient of the / y^— (yfi?^-terni in the classical action. 
One has to notice that this 9-order fine-tuning of the coefficient of the classical action of 
vacuum has confine this coefficient not to zero, but to some non-zero value c/(47r)^ ?« 0, so 
that the sum of two would be very close to zero. In our case of stable inflation we can take 
the coefficient of the / y/—gR^-teTm to be as small as ai, that does not spoil inflation and is 
in agreement with the renormalization group ^6 . 
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Figure captions 



Figure 1: Behaviour of tlie spectrum of CMB anisotropy for gravitational waves with k = 

-0.01. 

Figure 2: Behaviour of the spectrum of CMB anisotropy for gravitational waves with 
k = -0.2. 
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